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In general the kernel of QCD's gap equation possesses a domain of analyticity upon which the 
equation's solution at nonzero chemical potential is simply obtained from the in-vacuum result 
through analytic continuation. On this domain the single-quark number- and scalar-density distri- 
bution functions are /i- independent. This is illustrated via two models for the gap equation's kernel. 
The models are alike in concentrating support in the infrared. They differ in the form of the vertex 
but qualitatively the results are largely insensitive to the Ansatz. In vacuum both models realise 
chiral symmetry in the Nambu-Goldstone mode and in the chiral limit, with increasing chemical 
potential, exhibit a first-order chiral symmetry restoring transition at /i « M{0), where M{p^) is 
the dressed-quark mass function. There is evidence to suggest that any associated deconfinement 
transition is coincident and also of first-order. 

PACS numbers: 25.75.Nq 11.30.Rd 24.85. 4-p 12.38.Lg 



I. INTRODUCTION 

QCD possesses a C/b(1) symmetry, which is associated 
with a conserved charge; namely, baryon number, Qb- 
The finite baryon density theory is defined through the 
inclusion of a chemical potential, fj,, that is conjugate to 
Qb ■ The baryon number density is given by 



p(m) 






(1) 



where P{iJ.) is the thermodynamic pressure. Naturally, in 
the confined domain the only contribution to the pressure 
from an excess of quarks over antiquarks is that generated 
by the presence of colour-singlet baryon bound-states. 
That changes upon deconfinement. 

Confinement is a key emergent phenomenon in QCD. 
Another is dynamical chiral symmetry breaking (DCSB), 
which is responsible, amongst many other things, for the 
large mass splitting between parity partners in the spec- 
trum of light-quark hadrons, even though the relevant 
current-quark masses are small. Neither of these phe- 
nomena is apparent in QCD's Lagrangian and yet they 
play a dominant role in determining the observable char- 
acteristics of real-world QCD. 

In-medium; namely, for nonzero temperature and/or 
chemical potential, QCD will exhibit deconfinement and 
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chiral symmetry restoration. It is necessary to deter- 
mine at which values of the intensive thermodynamical 
parameters these transitions take place, whether they are 
coincident, and to elucidate the response of the order pa- 
rameters and also hadron properties to changes in the 
intensive parameters. That information is important, for 
example, in constructing the equation of state and un- 
derstanding the nature of matter in the core of dense 
astrophysical objects jl|. 

One might look to numerical simulations of lattice- 
regularised QCD for this information. Such studies sug- 
gest that deconfinement and chiral symmetry restoration 
are coincident at T ^ j2j. However, a satisfactory algo- 
rithm for nonzero chemical potential is lacking [3| . More- 
over, lattice-QCD has hitherto yielded little information 
on the in- medium behaviour of hadron properties [4] . 

A natural starting point for an exploration of DCSB is 
QCD's gap equation. There are numerous studies of this 
basic Dyson-Schwinger equation in-vacuum [5|, |6[ and in- 
medium [7]. It is noteworthy that extant self-consistent 
studies of concrete models of QCD, which exhibit both 
confinement and DCSB in vacuum, possess coincident de- 
confinement and chiral symmetry restoration transitions; 
e.g., Refs. |3, @, [l3, llllj [l2l] ■ This result appears to follow 
from the crucial role played by the in-medium evolution 
of the dressed-quark self-energy in both transitions. 

Confinement herein is generally understood to be ex- 
pressed through a violation of reflection positivity by 
coloured Schwinger functions [l3[ . This is known to be a 
sufficient but not necessary condition. 

On the other hand, there is no ambiguity about DCSB. 
It is signalled by the appearance of a gap equation solu- 



tion in which the Dirac-scalar piece is nonzero in the 
chiral hmit. This effect owes primarily to a dense cloud 
of gluons that clothes a low-momentum quark [Ij, [T^ . 
DCSB is the single most important mass generating 
mechanism for light-quark hadrons; e.g., one can iden- 
tify it as being responsible for roughly 98% of a proton's 
mass. A system in which the gap equation's DCSB solu- 
tion is the favoured ground state is said to realise chiral 
symmetry in the Nambu-Goldstone mode. The antithesis 
is termed a Wigner-Weyl realisation of the symmetry. 



Our theme is /j, 7^ because, compared with nonzero 
temperature, this is the poorer understood domain. 
Models and truncations of QCD can therefore provide im- 
portant information. Simple models suggest that a first- 
order chiral symmetry restoring phase transition occurs 
at /i = /i^ ^ Mq , where Mq is a constituent-quark mass- 
scale. This happens, e.g., in Refs. [3, [ifl], which exhibit 
a coincident deconfinement transition. It is noteworthy 
that below ^} the structure of the ground state can alter 
gradually but its qualitative character is unchanged. 



Our study complements existing work toward under- 
standing /z ^ in QCD. In addition to those noted above, 
kindred studies of: chiral symmetry restoration are de- 
scribed in Refs. [m jlTI . [ig . [l9| ; th e response of hadron 
properties, in Refs. [23, [U, [22, [23]; the possible realisa- 
tion of colour superconductivity, in Refs. [2l, [2E US, 123 ; 
and quark stars, in Ref. [2l,[2y|- One novelty of our study 
is a consideration of the effect of dressing the quark-gluon 
vertex. 



II. FORMAL STRUCTURE AND ANALYSIS 

Herein we elucidate the influence of /i on the nature of 
the solutions to QCD's gap equation, which reads"'^ 

Sip;n)-^ = Z2i^i■p+i-f4{p4 + ^^i)+m^'"') + ^ip; ^i) , (2) 

with the renormalised self energy expressed as 

.A 
S(p;/i) = Zi g'^{^i)Dp„{p~q;^) 

Jq 

XY7p5(<z;M)C(9,p;/i), (3) 

where J represents a translationally invariant regu- 
larisation of the integral, with A the regularisation 
mass-scale, Dp„{k;^i) is the dressed-gluon propagator, 
r°(g,p; /i) is the dressed-quark-gluon vertex, and tti^™ is 
the A-dependent current-quark bare mass. The quark- 
gluon-vertex and quark wave function renormalisation 
constants, Zi_2(C^j^^); depend on the renormalisation 
point, Q, the regularisation mass-scale and the gauge pa- 
rameter. A nonzero chemical potential introduces no di- 
vergences in addition to those present in the /Lt = the- 
ory. Hence, the renormalisation constants determined at 
/^ = are completely sufficient in-medium. 

At /i = the gap equation's solution has the form 



~s{py 



i-i-pA{p\e)+B{p\e)- 



(4) 



and the mass function M{p'^) = B{p'^,C,'^)/ A{p'^,C,'^) is 
renormalisation point independent. (Hereafter, when 
emphasis is useful or necessary, an overscore denotes a 
quantity calculated in-vacuum; viz., with /i = 0.) The 
propagator is obtained from Eq. Q augmented by the 
condition 



S{P) \2^^2 



i-fp + m{C) : 



(5) 



where 



This article is organised as follows. In Sec.[TT| we de- 
scribe the n ^ gap equation and its solution. This 
enables us to explain the notion of a domain of analytic- 
ity in /i, upon which the in-medium Schwinger functions 
can be obtained from their vacuum values via straightfor- 
ward analytic continuation, and the number- and scalar- 
densities, e.g., are /^-independent. In Sec. lIIII we intro- 
duce our models for the gap equation's kernel. The pri- 
mary model can be described as lying within the class 
of rainbow truncations of the gap equation. It can be 
systematically improved. The second differs in that we 
employ a vertex Ansatz, which is interpreted as modelling 
effects of dressing. Section HVl is extensive. It reports nu- 
merical results obtained in the chiral limit and includes 
a brief discussion of confinement ex pres sed through the 
violation of refiection positivity [iJ, [30, [U, |3^ ■ We re- 
capitulate and conclude in Sec.lVl 



m(C) = [Z2(C^A^)/Z4(C^A^)]m^-(A) 



(6) 



is the renormalised mass, with Z4 the Lagrangian mass 
renormalisation constant. The chiral limit means 771 — 
0, where rh is the renormalisation-group-invariant light- 
quark current-mass 33]. 

Poincare covariance is lost at /i 7^ 0. In this case the 
gap equation's solution can assume the general form [Sj] 

S{p;fJ.y^ = i'i-pA{p'^,p-u,C^) 
+ ij4iP4 + ifi)C{p^,p-u,C) + B{p^,p-uX^), (7) 



^ In our Euclidean metric: {7p,7(t} 



74717273; a ■ b 
like => p2 < 0. 



Ei=i o-ibi 



ll 



"^i—lO^ibi', and Pp time- 



where we have written u = (0, i/i). Rotational invariance 
remains. The propagator itself can be written 

S{p;fi) = ij-paA{p'^,p-u,0+crB{p'^,p-u,C'^) 

+ ij4{p4 + ifj,)ac{p'^,p-u,C'^)- (8) 

Consider now the kernel of the gap equation; viz., 
the integrand in Eq. ([3]), and chemical potential in the 
neighbourhood oi fi — 0. So long as the kernel is an- 
alytic on a set of nonzero measure in the complex-gl 
plane that includes the nonnegative real-g| axis, then 
the dressed-quark propagator is obtained from that in- 
vacuum through a straightforward analytic continuation. 
As recognised in Ref. [1^ , this procedure remains valid so 
long as no singularity in the kernel is contained or moves 
within the contour 



7 = limi^,^oo Uf^i7i ) 

71 = {z{t) ^t,-R<t<R}, 

72 = {z{t) = R + it,0<t<fi}, 

73 = {z{t) = t + ifi,R>t>-R}, 

74 = {z{t) = -R + it,fi>t>0}. 



(9) 



If one evaluates the integral in Eq. ^ using a free 
gauge boson propagator, a bare vertex and a free fermion 
propagator with mass M , then the conditions above are 
satisfied for fi < M . Plainly, this is a statement that 
nothing occurs in perturbation theory unless the chemi- 
cal potential exceeds the Fermi energy. 

On the other hand, suppose that nonperturbatively 
3/ia such that for fi < ^a the analyticity condition is 
satisfied, then 



V/i < fia ■ S{p; n) = S{p + u) , 
from which it follows that 

A{p^,p-u) = C_{p^,p-u) 
V/i<Ma : <( = A{{p + uf), 

B{p^,p-u) = B{{p + uf). 



(10) 



(11) 



DSE models have been explored in which /ia = [^, llOl 
and fia 1^ but small [l^ . Variants of the Nambu-Jona- 
Lasinio model can be constructed in which fia > Mq, 
where Mq is the model's constituent-quark mass-scale 

The gap equation is one of a system of Dyson- 
Schwinger equations (DSEs) [a,[3,[l3|, from whose struc- 
ture one may infer that 3 /xj > 0, an infimum, such that 
V/i < iM^ every Schwinger function can simultaneously be 
obtained via analytic continuation of its in- vacuum form. 
On this domain, in addition to Eq. PTI| one has, e.g.. 



Dp^{k;fi) ^ Dp„{k) , 
r;(p,g;/i) = fl{p + u,q + u). 



(12) 
(13) 



Another two important values of the chemical poten- 
tial are those above which quarks are no longer confined 
and chiral symmetry is no longer realised in the Nambu- 
Goldstone mode. We denote them, respectively, by /i^ 



and /x^. There are models in which /i^ — fx^ 0, [To| . 
That might also be the case in QCD. We write 



fie = inf(/i^,/i^) 



(14) 



It is noteworthy that the analytic properties of the 
dressed-quark propagator, for example, change abruptly 
at fic- 

It is relevant to list a few additional results that obtain 
on the domain < fj, < fM^, 



We define the quark number density 

(fip 



n^ifi) = 2NcNfZ2 

/i(b1;M) 



(15) 



/i(|p1;Ai), (16) 



1 

in 



(27r)3 
dp4tr-D{-j4:)S{p;fi). (17) 



For fi < fu^ one can be confident that Eq. (fTO|) is satisfied 
and hence 

-ifiKfil: hM-fi) = {) ^ n,{fi)={)- (18) 

namely, the system exhibits no excess of quarks over an- 
tiquarks.^ Plainly, on this same domain there is no con- 
tribution to the thermodynamic pressure from an excess 
of quarks over antiquarks. 

The vacuum quark condensate can be expressed [3a | 



dPp 



- {qq)l{fi) = 2NcZ, I j^ h{\p\;fi) , (19) 

/2(b1;M) 



1 

An 



(27r)3 
dpitroSip;fi). (20) 



The reasoning used above leads to the conclusion that 

Vm<M?: /2(b1;M) = /2(b1) ^ {qqf^ifi)^W/l; 

.(21) 
i.e., on this domain the vacuum quark condensate is in- 
dependent of fi. 

It will be evident at this point that the simplest situa- 
tion to understand is that with fic < /ij. As we have just 
explained, in this case the system's ground state is insen- 
sitive to increasing fi until the phase boundary is encoun- 
tered. We observed above, however, that instances exist 
in which fia < fic- Such are possibly the most physically 
relevant cases, in which order parameters and observables 
may exhibit a qualitative change at fj, = fia and thereafter 
evolve continuously on the domain fj, G {fia,fJ-c)- This is 
illustrated by Ref. 9], in which fia — and nq{fi) = on 
fi G [0, fic) but the condensate grows with increasing fi. 



^ For a noninteracting single-species Fermi gas this situation per- 
tains for n < M, where M is the fermion mass. 



III. GAP EQUATION MODELS 

Significant effort continues to be expended on deter- 
mining ttie precise nature of tlie kernei of QCD's gap 
equation. A diaiogue between DSE studies and resuits 
from numerical simulations of lattice-regularised QCD is 
providing important information; e.g., Refs. [33, [33, [33, 
H, SH, H EE H SI • This body of work can be used to 
formulate reasonable Ansdtze for the dressed-gluon prop- 
agator and dressed-quark-gluon vertex in Eq. Q . That 
is how we proceed. The approach can be particularly 
fruitful in connection with QCD at /i > since little is 
rigorously known about that theory. 

In connection with employment of these Ansdtze it 
is worth remarking that the DSEs admit at least one 
nonperturbative symmetry preserving truncation scheme 
[39I l46l . [43], which has enabled theproof of numerous 
exact results P, IH [H, [H [Ml, [13, [H H HI. It 
also provides a starting point for the formulation of re- 
liable models that can be used to illustrate those re- 
sults and make predictions with readily quantifiable er- 
rors [6, 56, 51, 58, 59, M] • 

The Ansdtze are typically implemented by writing 



Zig''Dp^{p-q)Tl{q,p) 



A". 



g{{p ~ qf) D'^^^ip - q)-T,{q,p) , (22) 



wherein Dp^'^{£) is the Landau-gauge free gauge boson 
propagator, 0(£^) is a model effective-interaction and 
Ta-{q,p) is a vertex Ansatz. In one widely used approach 
to in- vacuum physics, g{f') is chosen such that the one- 
loop renormalisation group behaviour of QCD is pre- 
served and the vertex is written 



^a{q.P) ^la- 



(23) 



This is the basis for a renormalisation-grou p-im proved 
rainbow-ladder truncation of QCD's DSEs P. [Sg. WX. [SSl. 
[59, 60]. 

If one is interested solely in the gap equation and 
not concerned with a symmetry preserving treatment of 
bound-state properties, then Ansdtze for the vertex can 
be employed whose diagrammatic content is unknown. 
The class of such models that has hitherto been employed 
can be characterised by [6l| 



iT,{k, t) = i^Aik", e^) Ja + (k+ i)a 



^-/■{k + i)AA{k',e) + AB{k^e^) 



where 









+ F{f)]. 

F{e) 



B-p 



(24) 

(25) 
(26) 



with F = A,B; viz., the scalar functions in Eq. ^. This 
Ansatz satisfies the vector Ward-Takahashi identity and 



on < /i < /ij its in-medium form is obtained through 
Eq. ini). For /^ > fi^, Eq. ^ generalises to [H 

ir^ (fc, £; fi) = iJ^Aik, t, fi)j^ + i^c(k, i; fi)j^^ 

+ {k + i)a 



--f^-{k + £)AA{k,e;fi) 



'--f\\■Ck + i)Acik,i■,^l) + ABCk,l,^^) 



(27) 



where k = k + u, 
T,Fik,£;^i) 
ApCk,i;ii) 



F{e,ki;^i) + F{e,li;iji) , (28) 



F{e,k4;fi)-F{f,ir,n) 
fc2 - p 



U7 • u, 



(29) 
(30) 
with 



and we have defined 7 = 7 — U7 • u, j" 
m2 = 1. 

One has a Slavnov- Taylor identity for the quark-gluon 
vertex in QCD, not a Ward-Takahashi identity. Hence, 
Eqs. ([24| and ([27| are not necessarily an improvement 
over Eq. ([^5[) . A comparison between results obtained 
with the different Ansdtze is nevertheless useful in iden- 
tifying those which might be robust. It is noteworthy 
that our vertex Ansdtze do not exhibit kinematic singu- 
larities [631 ■ 

Herein we employ a simplified form of the 
rcnormalisatioii-group- imp roved effective interaction 
in Rcfs. [y, [56|, [57|, [58], [59|, [60|; viz., we retain only that 
piece which expresses the long-range behaviour (s = fc^): 



Qi^) 



'4Dse-^/^ 



(31) 



This is a finite width representation of the form intro- 
duced in Ref. |£^], which has been rendered as an inte- 
grable regularisation of 1/fc^ (63. Equation (|?T|) delivers 
an ultraviolet finite model gap equation. Hence, the reg- 
ularisation mass-scale can be removed to infinity and the 
renormalisation constants set equal to one. We specify 
the fJ. ^ Ansatz via Eq. p2p . In making this modest and 
practical simplification we neverthele ss p roceed beyond 
the fi^O studies of, e.g., Refs. [i,[ll[2a[2ll,[2l[2l,[2l]. 
Studies arc underway with the complete interaction of 
Rcfs. ^,M,.5Z,.58„.59„M]• 



IV. NUMERICAL RESULTS 

The active parameters in Eq. pip are D and lu. How- 
ever, they are not independent. In reconsidering a 
renormalisation-group-improved rainbow-ladder fit to a 
selection of ground state observables [531 , Ref. [53] noted 
that a change in D is compensated by an alteration of w. 
This feature has further been elucidated and exploited 
in Ref. Wf\. With the interaction specified by Eqs. 



TABLE I: Results obtained for selected quantities with the 
vertex and D parameter value indicated: ^(0), Af(0) are 
p = in-vacuum values of the scalar functions defined in 
connection with Eq. Q ; the vacuum quark condensate is de- 
fined in Eq. (|19|l : and M, T are discussed in connection with 
Eq. (|40|) . All calculated quantities in GeV except D, in GeV'^, 
and A(0), dimensionless. Unless otherwise noted, all calcula- 
tions reported herein were performed in the chiral limit. 



Vertex 


D 


A(0) 


M{Q) 


-m)y" 


A/i 


Ti 


M2 


Ta 


Eq.dlSl) 


1 


1.3 


0.40 


0.25 


0.53 


0.17 


0.82 


0.0 


Eq. p4l) 


1 
9 


1.1 


0.28 


0.26 


0.30 


0.0 


0.80 


0.20 



([23|) and ((3T|) . fitted in-vacuum low-energy observables 
are approximately constant along the trajectory^ 



ujD= (0.8 GeV)^ 



g ■ 



(32) 



Herein, we employ oj — 0.5 GeV. 

It is not possible to explore whether the behaviour ex- 
pressed in connection with Eq. ([32]) is also realised with 
the interaction specified by Eqs. ([22]), ^ and fSI]) be- 
cause Eq. ([24]) cannot yet be used in the implementation 
of a symmetry preserving DSE truncation. Hence, we 
keep u) = 0.5 GeV and choose D = 0.5 GeV^, a combi- 
nation that produces the chiral-limit vacuum quark con- 
densate in Row 3 of Table [H 



A. Rainbow truncation 

It is straightforward to solve the gap equation with the 
interaction specified by Eqs. (P^ . ([^5]) and (pij) . the reg- 
ularisation mass-scale removed to infinity and the renor- 
malisation constants set equal to one. In Table |l] we re- 
port in-vacuum values characteristic of the solution. 

Our first novel goal was to determine /ia in Eq. pop . 
In a system that supports dynamical chiral symmetry 
breaking that can be done by analysing the following 
configuration-space Schwingcr function [66|: with t > 0, 



A(r,M) 



(2^ 



ip-x+ipAT 



5{p)<7B{p]iA- (33) 



Our results for a range of values of /i are depicted in 
Fig.[T] NB. In the absence of DCSB; i.e., when chiral 
symmetry is realised in the Wigner mode, gb = in the 
chiral limit. In the Wigner mode 










(34) 



because there is no dynamically generated mass-scale 
that can raise the Fermi energy of the massless theory. 



1. Reflection positivity 

In order to explain Fig. [1] we must first digress. Recall 
that in a quantum field theory defined by a Euclidean 
measure [30|, the Osterwalder-Schrader axioms [3l|, [S^l 
are five conditions which any moment of this measure 
(n-point Schwinger function) must satisfy if it is to have 
an analytic continuation to Minkowski space and hence 
an association with observable quantities. One of these 
is "0S3," the axiom of reflection positivity, which is vi- 
olated if the Schwinger function's Fourier transform to 
configuration space is not positive definite. The space of 
observable asymptotic states is spanned by eigenvectors 
of the theory's infrared Hamiltonian and no Schwinger 
function that breaches 0S3 has a correspondent in this 
space. Consequently, the violation of 0S3 is a sufficient 
condition for confinement. This connection has long been 
of interest and is explained, e.g., in Sec. 2 of Ref. |15| . 

For a free fermion of mass M 



S^'^'^p) = 



1 



ij-p + M 



(35) 



and one finds 



A^'^'^iT, fi) = - e-^^^'^^f'^eiM - /i) , (36) 

which is positive definite for /i < M. Hence, 0S3 is satis- 
fied on this domain and, indeed, there is certainly a well- 
defined single-particle asymptotic state associated with 
this Schwinger function. At 11 — M , on the other hand, 
this single-particle Schwinger function vanishes. That 
happens because the Dirac mass gap has vanished and 
the Fermi energy is breached. At this point an insta- 
bility appears in the theory, one that is associated with 
real-particle production. Plainly, as we indicated above, 
/id = AI for a free fermion. 

Suppose instead that one encounters a theory in which 
the dressed-quark propagator assumes the form [63, [63, 

m 



Scc{p) = - 



1 



1 



i^ ■ p + z ij ■ p + z* 



M + iF , 



(37) 

a function with conjugate poles in the coniplex-p4 plane; 
viz., poles at ±{iM ± F) when p = Q. In the absence of 
a chemical potential, this yields 



Aec(T) = ie-^^-cos(Fr), 



(38) 



a function that possesses infinitely many, regularly 
spaced zeros. Hence, 0S3 is violated and thus the 
fermion described by this Schwinger function has no cor- 
respondent in the space of observable asymptotic states. 
When fJ, j^ one has 



The value of rug is typical of the mass-scale associated with non- 
perturbative gluon dynamics. 



Acc(t, p.) 



1 



-(M-M)r cos(Fr) 6i(A// - /i) 



(39) 



This is a single-particle Schwinger function with a mass 
gap that vanishes when the chemical potential exceeds 
^{z = M + iV) = M, which is the magnitude of the 
imaginary part of the pole's position in the complex-p4 
plane. 

It is conceivable that a theory might produce 



N 

Sn{p) = J2t 



1 



ij ■p + z^ i"f ■p + z* 



(40) 



where, without loss of generality, < 5R(zi) < ^{zi+i), 
i — 1,2, . . . , N — 1; namely, a single- fermion Schwinger 
function possessing a sequence of N real and/or complex- 
conjugate poles with associated residues r.;. If this theory 
preserves the ultraviolet behaviour of QCD, Eq. ([5]), then 



N 



N 



^r,(C) = l,^r,(C)5ft(z,0-m(C). (41) 



1=1 
In this case 



i=l 



N 



i=l 



AAr(T,A.) = ^^e-W--)-^)-cos(5(z,)T)e(3?(zO-M): 



(42) 

which at fixed r > is a piecewise continuous function 
of /i. Plainly, one can therefore identify 



lia = ^{zi) : 



(43) 



i.e., with the location of the first discontinuity. Equa- 
tion (|40)l characterises a class of single-fermion Schwinger 
functions that contains all the cases we've considered. 
Hence, Eq. (|^5)) is a general statement valid for all theo- 
ries in this class. 




t[GeV"'] 



— -12 



H=0.53GeV 

^=0.54GeV 

^=0.55GeV 



10 



20 



30 



40 



t [GeV" 



2. Domain of analyticity 

We now return to Fig.[TJ The curves coincide with 
a Schwinger function of the types discussed in Eqs. ((571) 
and (HDl) . The evolution with r is dominated by the pole 
for which ^{z) is the smallest, and our fit yields zi = 
Ml + iTi, with the values given in Table HI The quality 
of the fit is good. In this case. 



^^a 



'}" = 0.53 GeV. 



(44) 



The lower panel depicts results obtained with /i > /i^*^. 
They are dramatically different from those in the upper 
panel. It will be plain from Sec. lIV AD that the evolu- 
tion with r indicates the presence of another singularity, 
Z2 — M2 + 1^2 in Eq. (1^0)) . The fitted values are given in 
Table |T1 No imaginary part appears in this instance. 

In Fig. [21 we present the evolution with /j, of the scalar 
functions in the dressed-quark propagator evaluated at 
p = 0. Consistent with the preceding discussion; viz., 
Eq. Uni), A{p = 0;p) = C{p = 0;/i) for /z < Ma''- 
Nonetheless, all the plotted quantities evolve with /z. 



FIG. 1: Schwinger function in Eq. (|33|) computed from the 
gap equation defined via Eqs. (|22|l . (|23|l and (|3ip . Upper 
panel - from bottom to top, /i = 0, 0.4 GeV, 0.5 GeV. Solid 
curves: our numerical results. Dashed curves: fit obtained 
via Eq. (|39|l . with M, V given in Table U Lower panel - 
fj,a <f^ = 0.53, 0.54, 0.55 GeV. 



In particular, M(0;/i) is monotonically increasing on 
this domain, as it is in all models that preserve the 
momentum-dependence of dressed-quark self-energies [7[ ■ 
Figure [3] displays the behaviour of the chiral- limit vac- 
uum quark condensate, Eq. p^ . Despite the fact that 
the scalar functions in S(j); fi) evolve with ^, the conden- 



sate is constant on /u < /u^ 



3. Phase transition 

In any study relating to QCD at nonzero chemical po- 
tential it is crucial to determine the point at which the 
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FIG. 2: Evolution with /i of the dimensionless quantities 
A{p = 0; n), C{p = 0; /i), and M{p = 0; /x) measured in GeV, 
all computed from the gap equation defined via Eqs. ([22|l . ([23|) 
and (I3T1). 



FIG. 4; Difference in pressure, Eq. (|47|) . between the Nambu- 
Goldstone and Wigner solutions of the gap equation defined 
through Eqs. (gUl, ^ and ([31]) . It is noteworthy that B(0) = 
(0.16 GeV)^ since a typical phenomenological value for bag 
constant is ~ (0.15 GeV)''. 
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FIG. 3: Chemical potential dependence of the vacuum quark 
condensate computed from the solutions of the gap equation 
defined via Eqs. im, ^ and ((STJ- 



theory makes a transition between the Nambu-Goldstone 
and Wigner realisations of chiral symmetry. All reason- 
able models support both these phases and those con- 
sidered herein are not exceptions. An exposition of the 
nature of these phases is presented in Ref. [70| • 

In a consistent rainbow truncation the transition can 
be studied by considering the pressure owing to dressed- 
quarks obtained through the "steepest descent" approx- 



imation; namely, 

ps(A*)=TrLn[5-i]-iTrp^], (45) 

where S is the solution of the gap equation, E is the as- 
sociated self energy, and "Tr" and "Ln" are extensions of 
"tr" and "In" to matrix- valued functions. Equation ([^5]) 
is just the auxiliary field effective action |71j , which yields 
the free fermion pressure in the absence of interactions; 
i.e., when S = 0. Owing to Eq. (|12p . in this analysis we 
can neglect the gluon contribution to the pressure. 

A system's ground state is that configuration for which 
the pressure is a global maximum or, equivalently, the 
effective-action is a global minimum. In the steepest- 
descent approximation the pressure difference between 
the Nambu-Goldstonc and Wigner phases is given by 



6(/i) 



AN, 






(27r4) 



PS, 

In 






PlC^ + B^ 



-\P\' K 



-2/2 



PIC? 



w 



'Cw 



(46) 



(47) 



wherein the subscript W distinguishes the solutions ob- 
tained in the Wigner phase from those obtained in the 
Nambu-Goldstone phase. Our result for this chemical- 
potential-dependent bag constant [721 ] is depicted in 
Fig-Hi It is evident that the Nambu-Goldstone phase 
is favoured on 



/z</.r', A^r' = 0.38GeV<Ma' 



rb 



(48) 

A quantitatively similar value of fJ.^'^^ is obtained, e.g., 
in the studies of Refs. |^, 0, [il|. 
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ployed^ to show that <ja„{p) satisfies 083 on /z < ^a- 
This is an indication that a first-order deconfinenient 
transition is coincident with chiral symmetry restora- 
tion in the rainbow truncation model under consider- 
ation, which was also the case in the DSE studies of 
Refs. gilii. 

In the upper panel of Fig. [5] we plot the single- 
particle number-density distribution function introduced 
in Eq. (fTTl) . It is identically zero for ij, < fi^'^^ < fia- For 
/i > ^J'*^ it is evaluated with the gap equation's Wigner 
solution and, since fi}'^^ > /i^, the result is nonzero. 
(The number density is discontinuous at /i^'"^.) Evi- 
dently, on this domain the model supports an excess of 
quarks over antiquarks, a feature which is consistent with 
deconfinement because baryon bound states are not in- 
corporated. On the Wigner domain one has Aw ^ Cw- 
Hence, as apparent in the figure, the Fermi momentum 
is 




PF = ^^[Cw{0)|Awm~0■^^J^■ 



(49) 
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FIG. 5: Upper panel - Single-particle number-density distri- 
bution function, /i(|p|,/x) introduced in Eq. (tlT)) . computed 
at /x = 0.38, 0.42 from the Wigner solution of the gap equa- 
tion defined through Eqs. (I22j, (l23j and ((3T|. This function 
is identically zero for /i < /ij''''. Lower panel - Single-particle 
scalar- density distribution function, f2(\p\,n) introduced in 
Eq. pOf) . computed a.t ^ = 0.3 from the Nambu-Goldstone 
solution of the aforementioned gap equation. This function is 
/i-independent for fi < fi^^^ and identically zero otherwise. 



The last of the observations in Eq. (|48p makes plain 
that in the model under consideration in this subsection 
an analytic continuation of the /i = solutions of the 
gap equation is sufficient to completely determine the 
phase structure. Moreover, in the chiral limit the tran- 
sition is first-order because the vacuum quark conden- 
sate is zero in the chiral-limit Wigner phase and hence 
the condensate changes discontinuously at fi — n'^'^^. 
This, toOjis in agreement with earlier DSE studies; e.g., 
Refs. [i,[l3,[ll,l23■ 
It is notable that the derivative criteria explained in 
Sec. 2 of Ref. h^ and exploited in Ref. fl2| can be em- 



The pointwise behaviour of the curves in Fig. [5] is a non- 
perturbative property of the model's Wigner phase in 
rainbow truncation, not a numerical artefact. Increased 
numerical precision does little more than sharpen the 
Fermi surface. Comparison with Fig. [H] indicates features 
that are sensitive to the truncation. The impact of a 
more realistic interaction [^, \5Q, [53, [5^, [5^, [60| is cur- 
rently being studied. 

The lower panel of Fig. [5] shows the single- 
particle scalar-density distribution function introduced in 
Eq. ([20]) . It is independent of fi on the domain fi < n^'^^ 
but vanishes in the Wigner phase. So long as /i^ < fia it 
is generally true that /2 vanishes when /i is nonzero and 
vice versa. 

The complete expression for the pressure will involve 
a contribution from colour-singlet baryon bound-states. 
These are expressed in the gap equation via corrections 
to the rainbow truncation. If a vertex Ansatz is employed 
whose diagrammatic content is unknowable, as has often 
been done and we do subsequently, one cannot quantify 
the baryon contribution. The genuine study of a baryon 
matter phase will only become possible when this sit- 
uation is improved upon. One anticipates it will yield 
a domain of chemical potential upon which the baryon 
number density and vacuum quark condensate are simul- 
taneously nonzero. 



B. Dressed vertex 

We now turn to a consideration of the dressed-vertex 
gap equation defined by Eqs. (HH), (HH) and ^^. The in- 
vacuum solution is characterised by the quantities listed 



* It is far easier to use this criterion than to compute the analogue 
of Eq. II33I I obtained with cah/ ■ 
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FIG. 6: Schwinger function in Eq. (|33|l computed from the gap 
equation defined via Eqs. (|22p . (|24|l and pi|l . Solid curves - 
top to bottom, ^ = 0, 0.2 GeV, 0.25 GeV. Dotted curve - 
fi = 0.31 GeV. 



in Row 3 of Table [H In Fig. [5] we plot the Schwinger 
function defined in Eq. (j33[) . which is evidently in the 
class epitomised by Eq. (|^^ . As usual the evolution with 
T at small-/! is dominated by the pole for which 5R(z) in 
Eq. (|40p is the smallest. In this instance our fit yields 
zi = Ml + iFi, with the values given in Table [H The 
quality of the fit is good and the absence of an imaginary 
part is noteworthy. In this model, 






0.30 GeV. 



(50) 



Kindred to Fig.[l] the pointwisc evolution of the 
Schwinger function with r is markedly different for /i > 
/Iq'^. It is consistent with the presence of at least one pair 
of additional singularities in Eq. (|10|) , which are charac- 
terised by Z2 = M2 + iT2 with the fitted values listed in 
Table HI For /j, = the presence of these other singular- 
ities does not necessarily entail that this dressed-quark 
Schwinger function violates the axiom of reflection pos- 
itivity. Hence, 0S3 alone cannot be used to argue for 
quark confinement in this instance. 

The gap equation associated with Eqs. (|^ . ([M]) and 
(PIj) certainly admits solutions that exhibit chiral sym- 
metry in the Nambu-Goldstone and Wigner modes. How- 
ever, owing to the fact that the diagrammatic content of 
Eq. (j24p is unknowable, it is impossible to write the truly 
corresponding dressed-quark contribution to the thermo- 
dynamic pressure. 

How, then, is one to determine the critical chemical 
potential for chiral symmetry restoration, /x^^'^? An es- 
timate can be made by considering the chiral susceptibil- 
ities 



xrn{^'^)=^BP{Q■^l) 
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FIG. 7: Chiral susceptibilities defined in Eq. (|51l) and com- 
puted from the gap equation defined via Eqs. ([22}, (|24|l and 
(|3ip . Upper panel - Nambu-Goldstone solution; lower panel 
- Wigner solution. 



where m is the current-quark mass and also a source- 
term for (qg), and P = NG or W indicates whether the 
scalar self-energy was calculated in the Nambu-Goldstone 
or Wigner phase [TSj.^ 

In the neighbourhood of /x = the Nambu-Goldstone 
phase is realised, a fact signalled by X^'^{ii) > and 
X^ {^) < 0. The Nambu-Goldstone susceptibility is pos- 
itive because the B^'^{p; fi) solution of the gap equation 
characterises the true vacuum and hence must be stable 
against a variation in the source-term for a chiral conden- 
sate. X^-^'-^ will remain positive until a point jjl — fJ,^^, 



(51) 



This procedure can work because the Wigner solution of the gap 
equation is accessible even for nonzero current-quark mass iTOll . 
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FIG. 8: Chemical potential dependence of the chiral- limit vac- 
uum quark condensate computed from the gap equation de- 
fined via Eqs. (I22j, ^ and ((3T|- Recall that /ijj'^ = 0.30 GeV. 



whereat the system becomes unstable against variations 
in the chiral-condensate source-term and X^^{^) there- 
fore exhibits a pole. On the other hand, at small fi 
the Wigner susceptibility is negative because the Wigner 
phase is unstable against the formation of a quark con- 
densate. As /i is increased, X^ {(x) will remain negative 
until a point ji — ji^ whereat the Wigner vacuum be- 
comes metastable and X^ {^) exhibits a pole. Plainly, it 
must be that in any given theory 



^r < M? < M^ 



(52) 



Moreover, based on the results in extant model studies 
|35l uM ' ^^ anticipate that 



^^i 



A^f ^ + A^r 



(53) 
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FIG. 9: Upper panel ~ Single-fermion number-density distri- 
bution function, /i(|p|,/i, computed from the Wigner solu- 
tion of the gap equation defined via Eqs. (|22|) . (|24p and (|3Hl : 
lower panel - single-fermion scalar-density distribution func- 
tion computed from this gap equation. 



As an additional example, for the rainbow truncation 
described in Sect. lIVA^ this expression predicts /ij = 
(0.54 + 0.24)/2 = 0.39 GeV, which compares well with 
Eq. dMI). 

In Fig.[7]we plot the susceptibilities defined in Eq. (jSTj) 
and computed from the gap equation defined via 
Eqs. dHj), (1211) and ^. They exhibit the anticipated 
behaviour and enable us to infer 



Aij'''«0.3GeV 






Evidently, in the dressed-vertex model we cannot at 
present separate IJ,^^'^ from ^^ . Hereafter we will as- 
sume these points to be coincident. 

The evolution with chemical potential of the chiral- 
limit dressed-vertex vacuum quark condensate is depicted 
in Fig. [SI It is constant for /x < ^^^''^ and vanishes for 



M > Ai: 



Xbc 



The chiral symmetry restoring transition is 



therefore first-order. Employing the derivative criteria 
explained in Sec. 2 of Ref. [1^ and employed in Ref. p^ . 
we find in the dressed-vertex case, too, that (JAw (p) sat- 
isfies 083 on /.t < /ia. 

In the top panel of Fig. [5] we plot the chiral-limit 
dressed-vertex fermion number-density distribution func- 
tion. It is only nonzero in the Wigner phase; i.e., for 
(54) M > Mc '"^ ill Eq. (|M)) . Therefore on this domain the pres- 
sure receives a contribution from an excess of quarks over 
antiquarks. This is consistent with deconfincmcnt but 
more cannot reliably be said because the diagrammatic 
content of the vertex Ansatz is unknowable. 

In the lower panel of Fig.[n] we depict the chiral-limit 
fermion scalar-density distribution function. It is inde- 
pendent of chemical potential for /i < /ij ^'^ and is only 
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nonzero in the Nambu-Goldstone phase. 



V. SUMMARY AND CONCLUSIONS 

The gap equation is the natural starting point for con- 
tinuum studies of chiral symmetry restoration at nonzero 
chemical potential; viz., fJ. ^ 0. We have explored and il- 
lustrated its features using two models for the equation's 
kernel that are distinguished from each other by the form 
of the vertex Ansatz. It is important that qualitatively 
our results are largely insensitive to the Ansatz. 

In general the gap equation's kernel possesses a domain 
of analyticity in fi. This means that there exists a /Za such 
that for all /x < /i^ the gap equation's solution is simply 
obtained from the in- vacuum result through analytic con- 
tinuation. While examples exist in which /!„ = 0, that is 
not necessarily the case and hence this observation can 
be practically useful. For example, it guarantees that 
when calculated with the gap equation solution which 
expresses dynamical chiral symmetry breaking (DCSB), 
the single-quark number- and scalar-density distribution 
functions are /^-independent on /x < /Iq . This will also be 
true of other physical quantities. 

A striking signature of dynamical chiral symmetry 
breaking is a nonzero value of the in-vacuum dressed- 
quark mass function, M(p^), which entails a nonzero 
chiral-limit vacuum quark condensate. Both models 
possess this feature and realise chiral symmetry in the 
Nambu-Goldstone mode. Naturally, a nonzero chemical 



potential destabilises the quark condensate and, with in- 
creasing chemical potential, both models exhibit a first- 
order chiral symmetry restoring transition at /x « A/(0). 
We described evidence which suggests the existence of a 
coincident first-order deconfinement transition. 

The models we studied are novel examples of the case 
M(0) ~ /^^ < Ma- However, they are alike flawed in the 
ultraviolet behaviour of the kernels. This can be reme- 
died by incorporating a term in the effective interaction 
which ensures that the one-loop renormalisation group 
behaviour of QCD is respected. A kernel of this type has 
been widely employed in the study of hadron properties 
in-vacuum [^, [5y, [53, S, [S^, [SO] • An exploration of the 
/i 7^ properties of such a model is underway. 
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